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I. INTRODUCTION 

Strain engineering is an efficient mechanical approach 
to manipulating the physical properties in quasi two- 
dimensional nanostructures such as graphene, M 0 S 2 , 
black phosphorus, and others. A huge number of works 
have been performed to examine the effectiveness of 
strain in modulating the physical properties of these 2D 
materials, with the particularly well-known example of 
using strain to generate a finite electronic bandgap for 
graphene (for a review, see eg. Ref. [l). 

Mechanical strain has also been used to modify 
the physical properties in single-layer black phosphorus 
(SLBP)— 9 In particular, it has been shown in a number 
of previous works that mechanical strain is an effective 
means to tune the electronic bandgap in a wide range 
for SLBP. A large uniaxial strain in the direction nor¬ 
mal to the SLBP plane can even induce a semiconductor- 
metal transition— c— The in-plane uniaxial strains along 
the armchair and zigzag directions have also been used 
to modify the bandgap of SLBPjIA— while the rela¬ 
tive efficacy of uniaxial and biaxial strains have been 
comparatively studied for their effects on the electronic 
band structure for SLBP—— First-principles calcula¬ 
tions have shown that both biaxial and uniaxial strains 
rotate the preferred electrical conducting direction by 90 
degrees— The method of invariants has been applied to 
investigate the electronic band structure of SLBP with 
external fields including the strain field*? 1 

However, in nearly all of the above works, the bandgap 
changes have been obtained through strains applied ei¬ 
ther in the armchair or zigzag directions, or in the di¬ 
rection normal to the SLBP plane. This is reasonable, 
because these three directions are principal directions 
for the D’ 2 h symmetry of the puckered configuration of 
the SLBP. 2 - However, a very recent study has demon¬ 
strated that the maximum in-plane Young’s modulus for 


the SLBP is neither in the armchair direction, nor the 
zigzag direction. Instead, there exists a third principal 
direction with direction angle <f> = 0.2687T, along which 
the SLBP has the largest Young’s modulus value— Sim¬ 
ilarly, there is no guarantee that the most effective mod¬ 
ulation of the bandgap by applying strain occurs in the 
armchair or zigzag direction. 

Furthermore, in most existing studies, the mechanical 
strain that is applied to the SLBP has been limited to 
either uniaxial or biaxial strain. Hence, a natural ques¬ 
tion to ask and answer is what the optimal direction and 
type of mechanical strain is that results in the largest 
variations in the bandgap. A systematic analysis and 
understanding for the strain effect on the bandgap for a 
general strain type will be essential for practical strain- 
based manipulation of the electronic properties in SLBP. 
This comprises the focus of the present work. 

In this paper, using the tight-binding approximation 
(TBA) model, we derive an analytic formula for the 
strain dependence of the electronic bandgap in SLBP. 
We obtain an analytic expression for the direction of the 
applied strain, along which the strain will induce the 
strongest modulation in the bandgap of the SLBP. In 
particular, the effects from different strain types (tension, 
shear, and coupled tension and shear) are systematically 
compared. 

The present paper is organized as follows. In Sec.II, 
we present details regarding the structure of SLBP. The 
TBA model for SLBP is introduced in Sec.III. Sec.IV (A) 
is devoted to the derivation of a general analytic for¬ 
mula for the strain dependent bandgap, and the bandgap 
variations induced by different strain types are compared 
in Sec.IV (B). The paper ends with a brief summary in 
Sec.V. 
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FIG. 1: (Color online) SLBP structure. There are two prin¬ 
cipal directions, i.e., the armchair (blue arrows) and zigzag 
(red arrows) directions. Color is with respect to the atomic 
z-coordinate. 


II. STRUCTURE 

The atomic configuration of the SLBP is shown 
in Fig. [0 The structure parameters were measured 
experimentally.— The two in-plane lattice constants are 
a i = r 3 7 = 4.376 A and a 2 = T 2 a = 3.314 A, while the 
out-of-plane lattice constant is a 3 = 10.478 A. The ori¬ 
gin of the Cartesian coordinate system is located in the 
middle of r*i 2 . The x-axis is in the horizontal direction 
and the y-axis is in the vertical direction. The z-axis is 
in the direction normal to the SLBP plane. There are 
four inequivalent atoms in the unit cell a\ x d 2 of SLBP, 
which will be chosen as atoms 1, 2, 3, and 6 in this work. 
The coordinates of these atoms are r 4 = (— ua\, 0, —va 3 ), 
r *2 = (ua\,0,va 3 ), r 3 = (0.5ai — ua\, 0 . 5 a 2 , va 3 ), and 
r *6 = (—0.5ai + ttai, 0 . 5 a 2 , —ua 3 ). The two dimension¬ 
less parameters are u = 0.0806 and v = 0.1017. The 
bond lengths from the experiment are d 3 = r 2 3 = r 16 = 
2.2449 A and d 2 = ri 2 = 2.2340 A, and the two angles 
are 0 3 28 = 0 . 5357 T and 0 3 2 i = 0.5677T. 


III. TBA MODEL FOR SLBP 

We describe now the electronic band structure for 
SLBP obtained using a two orbital TBA model, which 
is derived from a recently proposed four orbital TBA 
model— Specifically, it was proposed that the electronic 
band structure of the SLBP can be treated by a four 
orbital TBA model^ with four hopping parameters be¬ 
tween atom pairs (2, 3), (2, 1), (2, 6), and (3, 6) in 
Fig. □ Among these four hopping parameters, it was 


shown that the electronic band structure in SLBP is de¬ 
termined mainly by the first two nearest-neighbor hop¬ 
ping parameters between atom pairs (2, 3) and (2, 1). As 
a consequence, we use these two leading hopping param¬ 
eters to describe the electronic band structure for SLBP 
in the present work. 

The two hopping parameters in this two orbital model 
are t\ between atoms 2 and 3, and £2 between atoms 2 
and 1 . For undeformed SLBP, the hopping parameter 
between atoms 2 and 8 (t 3 ) is the same as that between 
atoms 2 and 3. After the SLBP is deformed by the me¬ 
chanical strain, hopping parameters <1 and t 3 become 
different, so generally we have three hopping parameters 
in the following. 

Based on the two orbital TBA model, the electronic 
Hamiltonian for the SLBP can be written as, 

/ Eq t2 0 t 3 T tl&2 

TT _ t 2 Eq tl + t 3 &2 0 

0 tl + t 3 S 2 Eq t 2 5i 

V t 3 -(- 1\8 2 0 t 2 S* x Eq 

where <5i = e lfclQl and 8 2 = e* fe2 “ 2 are two phase factors, 
with k = k\bi + k 2 b 2 as the wave vector. Here bi are two 
reciprocal bases defined by bi • dj = 2nSij for i, j = 1 , 2. 
The atomic energy level, Eq, is set to 0 in the following 
calculation. 

The eigenvalue solution for the Hamiltonian in Eq. ED 
gives four electronic bands for SLBP, 

C 4 E 4 + C 2 E 2 + Cq = 0, (2) 

where the coefficients are, 

C 4 = 1; 

C *2 = —2 (ti + f 2 +1 3 + 2fif 3 cos A 2 ); 

Co = t 2 ~ 2t2(cosAi [2tif 3 + (t\ +f 3 ) cosA 2 ] 

-sin Ai sin A 2 (t\ - f 3 )} + (t\ + t 3 + 2fif 3 cos A 2 ) 2 , 

where Ai = 2Ttk\a\ and A 2 = 2 Trk 2 a, 2 - 

Fig.[2]shows the electronic band structure from Eq. 0 
for undeformed SLBP along high symmetric lines in the 
first Brillouin zone. The two hopping parameters are 
t\ = —0.797 eV and = 2.393 eV. We have used the 
subscript 0 to denote hopping parameters in undeformed 
SLBP. For undeformed SLBP, the hopping parameter be¬ 
tween atoms 2 and 8 (t 3 ) is the same as that between 
atoms 2 and 3, i.e., t 3 = t\ = —0.797 eV. These param¬ 
eters are obtained from the corresponding hopping pa¬ 
rameters in the four orbital model by scaling them with 
the same factor, so that the bandgap from the two or¬ 
bital TBA model agrees with that from the original four 
orbital model.— The band structure in Fig. [2]is very sim¬ 
ilar as that from the four orbital TBA model, and in 
particular the conductance band and the valence band 
are very close to the four orbital model. The electronic 
band structure around the T point is an even function of 
wave vector as restricted by the symmetry group of the 
SLBP— 
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FIG. 2: Electronic band structure for undeformed SLBP using 
the two orbital TBA model. The bandgap A E gap = 1.60 eV 
is reached at the T point. 

We focus on the bandgap modulated by mechanical 
strain in the SLBP. In the linear deformation regime, the 
direct bandgap locates at the T point. The wave vector 
k = 0 at P point, so the four electronic energy states are 

Ei = (ti + t 3 ) — t2\E 2 = — (fi + £ 3 ) — t2] 

E 3 = (t\ + t 3 ) + t 2 \ Ei = — {t\ + t 3 ) + t2- 

The energy gap is, 

Eg&p = E 3 — E 2 = 2 (t\ + t2 + £ 3 ) ■ (3) 

For undeformed SLBP, we find that the bandgap E gap = 
1.6 eV, which agrees well with the four orbital TBA 
model and other first-principles calculations ] 25 ’ 26 


IV. STRAIN EFFECT ON ELECTRONIC 
BANDGAP 


direction angle <f> is determined starting from the x-axis, 
and so the armchair direction is for <f> = 0, while the 
zigzag direction is for </> = We perform a coordinate 
transformation, by rotating the x-axis in Fig. [T] to the 
strain direction = e x cos <f> + e y sin <f>. The coordinates 
for a vector in this new coordinate system become 

( Xfj, \ / cos (j) sin (j) 0 \ / x \ 

V<t> = -sin</> cos (f) 0 ) ( ?/ ] , (4) 

H) V 0 0 1 )\z) 

where ( x , y , z) is the original coordinate for the vector, 
and the subscript <f> is to denote quantities in the new 
coordinate system. In the new coordinate system, the 
coordinates are deformed by an arbitrary linear mechan¬ 
ical strain as 

x e \ / l + e x 7 0 \ [ x <l>\ 

Ve = 7 1 + 0 i /0 , (5) 

Ze J \ 0 0 1 + e z J \ J 

where 7 is the shear component, while e x , e y , and e z 
are normal strains. The subscript e in the coordinate 
is to denote quantities after deformation. We have de¬ 
coupled the z component from the other two in-plane 
components, considering the quasi-two-dimensional na¬ 
ture of the SLBP structure. The mechanical strain is 
applied by deforming the SLBP structure directly, which 
ignores the Poisson effect and does not account for sub¬ 
sequent structural relaxation. This treatment results in 
an oversimplification of the strain state, particularly as 
compared to those that would result in experiments, such 
as substrate bending and stretching, or pressing via an 
STM tip. However, the error is reasonably small for lin¬ 
ear deformation regime, as shown by comparison with 
prior DFT simulations in the following discussions. 

In the linear deformation regime, the bond length r 
can be expanded as a function of all strain components, 
e x , e y , e z , and 7 as 


A. General formula for strain modulated bandgap 

We now consider the strain effect on the electronic 
bandgap of the SLBP. The electronic bands for SLBP 
are composed of s and p orbitals^ Moreover, the hop¬ 
ping parameter ( t) between s and p orbitals depends on 
the bond length (r) a a 27 ’ 28 t cx It has been assumed 
that the principal directions of the two neighboring Wan- 
nier orbitals keep their orientation along the bond vector 
of the two neighbor P atoms, such that the angular de¬ 
pendence does not play a role in the strain effect on the 
hopping parameter. Instead, the strain effect is realized 
through changing the bond length. Thus, the applied 
mechanical strain can affect electronic states (including 
the bandgap) through modifying the hopping parameters 
in the TBA model. 

We consider the deformation of SLBP under a general 
mechanical strain in the direction with angle <f>. The 


dr dr 

r = r 0 + -z— £x + 

u€.x oe v 


= r 0 


c 2 / 


dr dr 

+ ^— e z + w—7 
de z d 7 

a z e z + a s j, 


( 6 ) 


where we have introduced a as the strain-related geomet¬ 
rical coefficients. Recalling the relationship between the 
hopping parameter and the bond length, t ix we get 
the strain effect on the hopping parameter, 

/ 2 2 2 2 \ 

t = t 0 1 - — a x e x - a y e y - —a z e z - a s ~/ (7) 

V ro r 0 r 0 r 0 J 


According to Eq. Q, the key ingredient is to compute 
the strain-related geometrical coefficients a for each hop¬ 
ping parameter. For the strain e x , we get the following 
geometrical coefficients for each hopping parameter ti, 


x dr 23 , 1 2 

Oh = k=0 = - x 23<[> 

de x r 2 3 
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1 2 

= — [(0.5 — 2w) ai cos</> + 0 . 5 a 2 sin</>] ; 
d\ 


<9*28 , 


1 


“ 3 - 5^=° - 

1 


= -j- [(0.5 — 2 m ) ai cos</> — 0 . 5 a 2 sin 0 ] 2 ; 
d i 


x <9*21 | 1 2 l/o ,\2 

a 2 = -x —k=o =—*210 = — ( 2 uai cos 0 ) . 
Ok *21 «2 


Here, a x is the coefficient corresponding to the hopping 
parameter t\. 

For the strain e y , we obtain the following geometrical 
coefficients, 


v <9 * 23 1 

«» = _|„. 0 = 

- 5T [ - (0 ' 5 - 

v dr 2 s , 

= ——e„=0 = 


<9e„ 


*28 


1 2 
= — [(0.5 — 2rt) ai sine/) + 0 . 5 a 2 cost/)] : 

«i 

<9*21 1 2 l/o • 2 

—— = —j /210 = y ( 2 wai sin </> . 

aey r 2 i v d 2 


jj _ 


For the e z strain, we get the following geometrical co- 


efficients, 






< 9 r 23 


1 2 


a{ = 

de z 

k=o - 

*23 



< 9*28 


1 2 


«3 - 

de z 

k=o - 

- z 280 — 

*28 

0 ; 


< 9*21 


1 2 

1 

a 2 - 

de z 

k=o — 

2 210 — 
*21 V 

d-2 


We can derive similar expressions for the geometrical 
coefficients, a a , corresponding to shear strain, 


<*1 = —* 2302/230 
Ml 

2 

= — [(0.5 — 2w) ai cos ^) + 0 . 5 m 2 sinc/i] 

Ml 

x [— (0.5 — 2 u) mi sin </> + 0 . 5 a 2 cos <j )\; 

s <9*28 2 

= —7-17=0 = -*2802/280 

* 28 


<9y 


= — [(0.5 — 2m) mi cos</> — 0 . 5 m 2 sint/i] 

Ml 

x [— (0.5 — 2 u) Mi sin </> — 0 . 5 a ,2 cos (f >\; 


-4e x 

-4e z 


L l /„® I I C 2 Qf 2 

5r (o 1 +a s ) + ^ r 

^1 / 2 I z\ , ^2 a 2 
- {ai + a ,) + — 


-4 e 
— 47 


r +0 f o v 

>(a? + a|) + ^ 


y di 
t° 


t°rv s 

1 / s , „s\ , t 2 u! 2 

J -(a 1 +a,) + — 


1 2 

eo = 

2/230 
*23 V 



ei = 

2u) Mi sin </> + 0 . 5 a ,2 cos </>] 2 ; 


1 2 

Zl 00 a. 

e 2 = 


After some algebraic manipulation, we get the strain in¬ 
duced modification in the bandgap, 

AAgap = e 0 e z + (d - 2 e 2 ) (e x + e y ) - 2 e 2 e cos ( 2 </> + ip ), 

( 8 ) 

(9) 

( 10 ) 

= 0.411 eV. 

( 11 ) 


where the parameters eo, ei and e 2 are as follows 
+0 

2t\a\ 


di 

2 

df L 


= 3.507 eV; 


[(0.5-2„)«] 2 -|] + f ” (2 ““ 1>2 




We have introduced the following two quantities in the 
above derivation, 


tani/> 


' 




= y(^x-ey) +(2q) 2 . 


( 12 ) 

( 13 ) 


Eq. © shows the variation in the bandgap induced by 
a general strain applied in the direction with directional 
angle </>. As can be seen from Eq. ([8]). the variation in 
the bandgap depends on the strain angle </> with period 
7 r. For a given strain ratio, tan ip = f 2 e , the maximum 
(or minimum) strain effect can be achieved, if the strain 
is applied in the direction with angle </> satisfying 

cos (2 (f> + if>) = ±1, 
which gives the strain direction, 


4> .71 

~2‘ +J 2’ 


(14) 


^ r21 = ~~ ^ 2102/210 = -J- (— 2 um! cos (f>) (2ua± sin </>). 
2 


2 r 2 r~ d: 


Inserting these geometrical coefficients into Eq. 0 , 
and using Eq. ©, we obtain the analytic expression for 
the strain dependence of the electronic bandgap, 


rp _ Z7>0 _ 

-£Agap -C'gap — 


where j is an integer. This means that mechanical strain 
can introduce the largest (smallest) modulation of the 
bandgap if the strain is applied in the direction described 
by Eq. (fill) . In particular, we note that, to achieve the 
strongest strain effect on the bandgap, there is no guar¬ 
antee that the strain should be applied in the armchair 
or zigzag direction. Instead, the optimal strain direction 
is generally dependent on the type of the applied strain. 


B. Comparison between different strain types 

In the above, we have derived the bandgap variation 
induced by a general strain in Eq. ©. We have also 
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obtained the direction for a general strain in Eq. m, 
where the direction lies in the 2D plane. This direc¬ 
tion represents the most efficient strain direction, in that 
strain applied in this direction will generate the largest 
modulation of the bandgap. 

In this section, we will determine the most efficient 
direction for some common strain types in SLBP. We 
first note that e\ — 2&i > 0 , ei > 0 , and e = 

\J (fix — e y ) * 1 2 3 4 + (2 7 ) 2 > 0 in Eq. ([5]). It is obvious that 
strains e x and e y have similar effects on the bandgap, so 
we will discuss only one of them in some situations in the 
following. 

(1) For uniaxial strain in the z-direction, i.e., e x = e y = 
0 , 7 = 0 , and e z ^ 0 , we have 

AF ga p = Co C z . (15) 

We can see that the change of the bandgap is a linear 
function of the applied strain. This is consistent with 
previous first-principles calculations— r— 

(2) For in-plane uniaxial strain, i.e., e x ^ 0, e y = 0, 
e z = 0 and 7 = 0 , we have 

AFg ap — (ei 2 e 2 ) € x ‘2 c2^ x COS 2<p. (15) 

The most effective direction is determined by the condi¬ 
tion that both terms on the right side have the same sign, 
i.e., cos2 cj) = —1. This gives cf> = which is the zigzag 
direction in SLBP, and means that uniaxial strain can in¬ 
troduce the strongest effect on the bandgap if it is applied 
in the zigzag direction in SLBP. For this uniaxial strain 
in the zigzag direction, the bandgap is AF gap = e\t x . 
The coefficient e\ > 0, leading to an increase of the 
bandgap due to tensile strain, which is consistent with 
first-principles calculations— 

As another example, if we assume that the uniaxial 
strain e x is applied in the armchair direction (c i> = 0 ), 
then we have AF gap = (ei — 462 ) 63 ,, where the coefficient 
(ei — 4 e 2 ) < ep. This means that, to induce the same 
bandgap variation, a larger strain magnitude is needed if 
the uniaxial strain is applied in the armchair direction. 

(3) For in-plane biaxial strain, i.e., e x = e y = e, e z = 0 
and 7 = 0 , we find, 

AF gap = (ep — 2 e 2 ) ( e x + e y ) = 2 (ei — 2 e 2 ) e. (17) 

There is no preferred strain direction for biaxial strain, 
which is consistent with the intrinsically isotropic nature 
of biaxial strain. 

(4) For a general in-plane strain with e x ^ e y , e z = 0 
and 7 = 0 , we find 

A-Fgap (Cl ^ 62 ) ( 63 , T 6 y) 262 ( 63 : 6 y) COS 2(j). 

(18) 

The most efficient strain direction depends on the sign of 
AFg a p. More specifically, it requires both terms on the 
right side to have the same sign as AF gap . 

For AFgap > 0, an effective strain application should 
require ea,+ey > 0 according to the first term on the right 



FIG. 3: (Color online) The most effective approach to en¬ 
larging the bandgap by a combination of normal strain e x 
and shear strain 7 . The direction angle for the strain is 
(f> = — ^ + (2j + 1 )-|, with tan ip = 


side. From the second term, we have —( 63 , — e y ) cos 2(f> > 
0 ; i.e., we should have e x < e y for ( f > = 0 or 63, > e y 
for (j) = This indicates that the tensile strain should 
be applied in the two principal directions (armchair and 
zigzag) of SLBP, so that the bandgap can be enlarged 
most effectively. Furthermore, for maximum bandgap 
increase, the tensile strain should be larger in the zigzag 
direction than the armchair direction. 

For AFg ap < 0, the most effective strain application 
for bandgap reduction should require e x + e y < 0 accord¬ 
ing to the first term on the right side. From the second 
term, we have — (e x — e y ) cos 2 </> < 0 ; i.e., we should have 
63, > 6y for < f > = 0 or 63, < e y for ( f > = ■|. This indi¬ 
cates that the axial strain should be applied in the two 
principal directions (armchair and zigzag) of SLBP, so 
that the bandgap can be reduced most effectively. Fur¬ 
thermore, the compressive strain should be larger in the 
armchair direction than the zigzag direction to reduce 
the bandgap. Considering that the strain is compressive 
in this situation, we have larger strain magnitude in the 
zigzag direction than the armchair direction. 

As a result, for both AF gap > 0 and AF gap < 0, 
strains should be applied in the two principal directions 
(armchair and zigzag) of SLBP. This result is consistent 
with recent first-principles calculations— Furthermore, 
the strain magnitude in the zigzag direction should be 
larger than the strain magnitude in the armchair direc¬ 
tion to achieve the largest bandgap change. 

(5) For pure shear strain, i.e., e x = e y = e z = 0 and 
7 ^ 0 , we find 

AFg a p = 4 e 2 7 sin 2 <)>. (19) 

It is important to point out that the most effective direc¬ 
tion for the shear strain is determined by sin 2 </> = ± 1 , 
which gives <f> = ± \, which illustrates that the most effec¬ 
tive direction for pure shear is not in either the armchair 
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TABLE I: Summary for the strain dependent bandgap variation. The last line lists the most effective direction for each strain 
type, along which the maximum bandgap variation can be achieved. 


strain type 

uniaxial strain 

uniaxial strain 

biaxial strain 

general strain 

shear 

uniaxial strain and shear 

definition 

6 X — £y = 0 

/ 0,7 = 0 

€x / 0? €y = 0 

e z = 0,7 = 0 

€-x €y C 

e* = 0,7 = 0 

€x / e y 
e z = 0,7 = 0 

€ x — €y — 0 
= 0,7/0 

€x / 0, 6y = 0 
=0,7/0 

Ai?g a p 

eoe z 

(ei - 2e 2 ) t x 
—2e 2 t x cos 2(f> 

2(ei — 2e 2 ) e 

(ei — 2e2) (e x + € y ) 
— 2 e 2 \e x — e y \ cos20 

2e 2 7sin 2 <f> 

(ei - 2e 2 ) tx 

—2e 2 yjtl + 7 2 cos ( 2(p + ip) 

0 max 

AE gap > 0 

N.A. 

zigzag, tf>= f 

arbitrary 

0 = §, e x > ty > 0 

-e- 

II 

H- 

d)=-± + (2j + l)% 

AEgap < 0 

(j) = 0, ty < tx < 0 

+ 

~ I" 

1 

II 

-e- 



FIG. 4: (Color online) The most effective approach to de¬ 
creasing the bandgap by a combination of normal strain e x 
and shear strain 7. 

or zigzag directions of SLBP. Instead, a pure shear strain 
should be applied in the direction with (j> = ±^, so that 
it can introduce the strongest effect on the bandgap for 
the SLBP. 

( 6 ) For strain with e y = e z = 0, e x ^ 0, and 7 yf 0, 
we simultaneously apply the uniaxial strain e x and the 
shear strain 7 to modulate the bandgap of SLBP. In this 
situation, we have, 

AE gap = (ei - 2e 2 ) e x - 2e 2 ecos(2</> + if>) . (20) 

To enlarge the bandgap, i.e., AE gap > 0, it can be seen 
from Eq. (EU1) that the most effective direction for apply¬ 
ing strain is to ensure cos (2 <f> + if)) = — 1. This deter¬ 
mines the angle for the strain direction, 

<t> = -| + (2j + l)|, (21) 

where j is an integer. Furthermore, e x and 7 are related 
to each other as, 

AE gap = (ei - 2e 2 ) e x + 2e 2 \/+ (2y) 2 . (22) 

Fig. [3] shows this relation between e x and 7 for different 
AE gap . Each curve in the figure indicates the most ef¬ 
fective way to generate the corresponding change in the 


bandgap. It is clear that e x < 0 is not a good choice, 
because it requires larger shear strain 7 . Hence, for 
AE gap > 0, the most effective way is to apply a strain 
with e x > 0 , along with an appropriate, non-zero choice 
of shear strain 7 . If larger e x is applied, then the re¬ 
quired shear component 7 is smaller. We note again that 
the strain direction (</>) is determined by the actual ap¬ 
plied strain e x and 7 , because of the relationship between 
cf> and if) in Eq. (I2T1) . and because tan^ = -A. 

Similarly, to reduce the bandgap, i.e. AE gap < 0, the 
most effective direction for applying strain is to ensure 
cos (2(f> + if>) = 1. This determines the angle for the strain 
direction, 

0=-| + 2jx|, (23) 

where j is an integer. Furthermore, the strains e x and 7 
are determined by the following relation, 

AE ga p = (ei - 2e 2 ) e x - 2 e 2 \j e 2 + (27 ) 2 . (24) 

Fig. [4] shows this relation between e x and 7 for different 
AE gap . 

The above discussions on different strain types are 
summarized in Tab. [T] From the third line in the table, 
uniaxial strain in the direction normal to the SLBP plane 
is the most effective strain type to modify the bandgap. 
In other words, to generate the same bandgap variation, 
this strain type requires the smallest strain magnitude 
among all strain types that have been discussed, because 
it has the largest pre-coefficient magnitude, |eo|. How¬ 
ever, the ability to apply different strain types, and com¬ 
binations of strain types, is dependent on the experimen¬ 
tal technique that is utilized. Thus, we expect that Tab.U 
can serve as a guideline for experimentalists to choose the 
most appropriate strain type to manipulate the bandgap. 

We note that all discussion in this work have been 
based on the TBA model, which does not consider the 
structural relaxation and orbital hybridization effects. 
This approximation is suitable for linear deformation 
regime, but not for nonlinear deformation with larger 
strains where structural relaxation and orbital hybridiza¬ 
tion occur. Because of this, the TBA model is not able to 
predict certain phemonena, such as the recently reported 
direct to indirect transition in the bandgap of SLBPjii 
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V. CONCLUSION 

In conclusion, we have developed an analytic model 
based on the tight binding approximation to elucidate 
strain effects on the electronic bandgap in single layer 
black phosphorus. We have demonstrated that the di¬ 
rection along which the mechanical strain is applied is 
critical to achieving the maximum modulation of the 
bandgap. More specifically, we have performed a detailed 
comparison between the effects from different strain 


types, and for each strain type, we present predictions 
for the most efficient direction for the mechanical strain 
as summarized in Tab. |U 
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